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Abstract 

This paper addresses analytical aspects of deterministic, continuous-time dynamical 
systems defined on networks. The goal is to model and analyze certain phenomena 
which must be framed beyond the context of networked dynamical systems, understood 
as a set of interdependent dynamical systems defined on the nodes of a (possibly 
evolving) graph. In order to advocate for a more flexible approach to the study of 
network dynamics, we tackle some qualitative problems which do not fall in this working 
scenario. First, we address a stability problem on a network of heterogeneous agents, 
some of which are of dynamic nature while others just impose restrictions on the 
system behavior. Our second context assumes that the network is clustered, and we 
address a two-level stability problem involving the dynamics of both individual agents 
and groups. The aforementioned problems exhibit lines of non-isolated equilibria, 
and the analysis implicitly assumes a positiveness condition on the edge weights; the 
removal of this restriction complicates matters, and our third problem concerns a graph- 
theoretic characterization of the equilibrium set in the dynamics of certain networks 
with positive and negative weights, the results applying in particular to signed graphs. 

Our approach combines graph theory and dynamical systems theory, but also uses 
specific tools coming from linear algebra, including e.g. Gersgorin discs or Maxwell- 
type determinantal expansions. The results are of application in social, economic, and 
flow networks, among others, and are also aimed at motivating further research. 

Keywords: graph; digraph; signed graph; hypergraph; network; dynamical system; stabil¬ 
ity; eqnilibrinm; differential-algebraic eqnation; flow; social network; economic network. 

AMS subject classification: 05C21, 05C22, 05C50, 05C82, 34D20, 90B10, 91D30, 94C15. 

1 Introduction 

The interactions between dynamical systems and networks have attracted the interest of 
many researches in the last decade. Needless to say, the stndy of dynamics taking place 
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on graph-theoretic structures has a long tradition in different application fields, including 
nonlinear circuits, water and gas networks, neural nets, fiows, coupled dynamical systems 
and, in particular, coupled oscillators and synchronization, etc., to name but a few. In the 
framework of what is now known as network science, dynamics is one of the key research 
directions. A great deal of research within this interdisciplinary field is directed to the 
dynamics of state variables defined on the nodes of a network; the topics involved include 
networks of coupled dynamical systems |2Hl |39l 091 [52] , consensus and agreement protocols 
0Ol 03] (50] [51] , controllability and related aspects [Ml [35l 03] , epidemics [161 EH] , diffusion of 
information, innovations, etc. [HUE], dynamics of social networks [131 ESI [36] or evolutionary 
dynamics [21 [32], among others. 

Much attention has also been paid to the study of network growth [21 El 011 02], percola¬ 
tion [38] and, more generally, to the analysis of the dynamic changes of the network itself (cf. 
[21EH1SZ]). These dynamic changes of the network structure have been studied in the last 
decades both in a deterministic and in a probabilistic context. The original ideas regarding 
these evolving networks [2] can be traced back to the seminal work of Erdos and Renyi on the 
evolution of random graphs [nnn]. Find more references in [21 EHl EH]. In a deterministic, 
continuous-time context, the changes in the structure of a network modelled as a weighted 
simple graph can been naturally described as a dynamical system in which the state variables 
are the graph weights or, equivalently, the adjacency matrix, in the understanding that a 
vanishing weight Wij models the absence of an edge between nodes i and j. This dynamical 
system can be defined as a one-parameter group of transformations of a given state space, 
or as a system of differential equations. Find details in 07] . 

So-called coevolutionary processes involve simultaneously dynamics taking place on a 
network and changes in the network structure; cf. [191 EHl ESI [2H1ES] and references therein. 
Indeed, it is natural to consider networks in which dynamical changes in the nodes coexist 
with the evolution of the network itself. With different names, these coevolutionary processes 
have been addressed in several application fields, including neural network training [5l [2H 
06], social network analysis [25] [36], economic models [2H], etc. See also [191 EH] and the 
bibliography therein. 

In spite of the large amount of literature addressing dynamical processes on networks, it 
seems that in a continuous-time, deterministic context, general approaches to the formulation 
and analysis of network dynamics are somehow too restrictive, being focused on one of 
just three paradigms: a) systems in which the dynamics takes places only on the network 
nodes; b) dynamics of the network itself; and c) coevolution processes combining both. The 
references cited above support this point of view. Certainly, in specific fields, dynamical 
systems on networks arise beyond these categories; flow networks or nonlinear electrical 
circuits, which systematically exhibit state variables not supported on the network nodes, 
are examples of this. But a general approach, studying abstract dynamics on networks 
beyond the three settings indicated above, seems to be lacking in the literature. 

This paper is intended to be a step in this direction, addressing certain dynamical prob¬ 
lems which cannot be framed in the working scenarios mentioned above. First, models for 
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network dynamics formulated in terms of a set of dynamical systems defined on the nodes 
of a graph somehow assume a homogeneous nature on all agents. By contrast, we address in 
Section |2] a stability problem involving a potential-driven flow on a network of heterogeneous 
agents, some of which are of dynamic nature while others just impose restrictions on the 
system behavior; this naturally leads to a differential-algebraic formalism. In Section |3] we 
will examine clustered networks, addressing a two-level stability problem which involves the 
dynamics of both individual agents and groups. Lines of non-isolated equilibria will naturally 
arise in the two problems just mentioned; the one-dimensional nature of the equilibrium set 
is implicitly related to a positiveness condition on the edge weights. In (possibly nonlinear) 
flow problems on networks which may include negative weights, the structure of the equilib¬ 
rium set becomes more involved: in Section H] we will examine such structure without the 
aforementioned positiveness assumption, applying in particular the results to flow networks 
on signed graphs. Finally, Section |5] compiles some concluding remarks. 

2 Heterogeneous agents and constrained dynamics 

Along the terms presented in the Introduction, the “dynamics on nodes” scenario for mod¬ 
elling continuous-time, deterministic network dynamics assumes that the dynamical process 
can be described by an explicit ODE such as 

x' = f{x), (1) 

where x is the state vector, which is defined by a set of variables supported on the network 
nodes. The topology of the network is assumed to be comprised on the vector held / and, for 
the sake of simplicity, we are ignoring control variables, explicit time dependences associated 
with system inputs, etc. A very common perspective assumes that the network edges model a 
coupling between a set of (otherwise independent) dynamical systems dehned on the network 
nodes (see e.g. [28l [39l SHI |52]), in a way such that only the states of adjacent systems are 
visible by each node. This means that the differential equations above read as 

x'i = fi{xi,Xj„...,Xj^^.^), ( 2 ) 

for the Ath node (i = 1,.. .n); here x* is the set of state variables dehned on node i, and 
ji,... ,jd{i) are the nodes adjacent to i. Note that for a fully connected (complete) network, 
(1^ makes no diherence to (II]) (this is the case for instance in the Kuramoto model; see e.g. 

H). 

Needless to say, the setting described above has a very broad application scope, as the 
literature discussed in Section [1] shows. However, in other cases this may be too restrictive 
a framework to model certain dynamical processes on networks. In particular, systems ([I|) 
and (|2|) do not capture the eventual existence of constraints restricting the whole dynamics 
of the network. These constrained dynamical systems are better framed in the theory of 
diherential-algebraic equations (DAEs; cf. [T^ 130] IM ] IT5]). 
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In this context, we focus in this section on networks in which two different sets of nodes 
(agents) coexist; the behavior of some of them are defined by explicit dynamical systems, 
whereas the others do not involve dynamics explicitly but are governed, instead, by algebraic 
(non-differential) restrictions. Letting y stand for the dynamic variables, which model the 
states of agents with explicit dynamics, and 2 ; for the algebraic variables, which correspond 
to agents without explicit dynamics, these systems can be written as a so-called semiexplicit 
DAE 


y' = f{y,z) (3a) 

0 = g{y,z), (3b) 

where / and g capture the topology of the network; an example can be found in (fTTj) below. 
2.1 Potential-driven flows and the graph Laplacian 

To fix ideas, consider a simple connected network in which a given resource or commodity 
flows among a set of n agents located at the network nodes; here “simple” means that 
the network has neither parallel edges nor self-loops. Let us first assume that all agents 
collect this commodity, and denote the amount of the collected resource at node i by x,. 
The flowrate (or flow) of this resource at a given edge j, connecting nodes i and k, will 
be denoted by Uj. Every edge is endowed with a reference direction; if the j-th edge is 
directed say from i to k, Uj > 0 (resp. Uj < 0) means that the resource flows from i to k 
(resp. from k to i). This drives the problem to the context of directed graphs, although it is 
worth emphasizing that all the results are independent of the choice of directions. Denote 
X = (xi,..., x„), u = ( mi , ..., Um)i where m is the number of edges. 

With these reference directions, the entries of the incidence matrix A = {aij) G of 

the resulting digraph read as 

{ 1 if edge j leaves node i 

— 1 if edge j enters node i (4) 

0 if edge j is not incident with node i. 

The continuity equations at the nodes can be then written as 

x' = —Am, (5) 

since the increase x' of the stored resource at node i equals the net flow entering the node, 
that is, —AiU with the sign convention above and Aj standing for the i-th row of A. 

System (jS]) is underdetermined. The way in which the flowrates u interact with the node 
state variables x may be defined from different perspectives; for instance, in a game-theoretic 
setting the agents’ strategies would define the flows, or in the framework of control theory u 
might be designed as to achieve a given goal. Here we will simply assume that u is designed 
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in order to get a fair (equal) distribution of the commodity among all agents. This can be 
achieved by setting the following potential-driven dehnition of the flowrates: 

u = A^x, (6) 

which is just a redistribution law in which the flow from i to k equals Xi — Xk- It is a trivial 
matter to recast (j5])-([6]) as 

x' = —Lx (7) 

where L = AA^ is the graph Laplacian matrix. Note that ([7]) has the form depicted in 
([2]), as a consequence of the identity L = —C -|- D, where C is the adjacency matrix of 
the graph (without directions) and D is a diagonal matrix whose Tth entry is the degree 
of node i. The dynamics of ([7]) has been analyzed in the context of so-called consensus 
protocols [iQl HSl [501E] , and the state variables x may be checked to converge to a common 
value which is the arithmetic mean of the initial values Xi(0),... ,x„(0). This means that 
the resource is indeed redistributed among all the agents in a way such that all of them 
asymptotically store the same quantity. 

2.2 Heterogeneous agents 

We now drive our attention to a different problem; given again a simple connected network, 
assume that only some agents, say type-1 agents (which w.l.o.g. are assumed to be those 
from 1 to r, with 1 < r < n) accumulate the aforementioned resource, again in an amount 
Xi for the Tth agent. The continuity equations for type-1 agents read as above, 

x[ = -AiU, i = 1,.. .,r. (8) 

The other agents, say type-2 ones (those numbered from r -|- 1 to n) are different and do 
not collect the commodity; instead, they are however intermediaries among type-1 agents. 
Therefore, the continuity equations at type-2 agents are 

0 = AiU, i = r + 1,... ,n. (9) 

Finally, we assume that the flowrates are still given by ([6]); note that by means of the 
variables x* type-2 agents now simply set up a reference value for the flowrates u. Denote 
y = (xi,..., Xr), z = (xr+i,..., x„) (and, for later use, yi = Xj, Zi = Xj+r) and split A as 



where Ai (resp. A 2 ) is the submatrix of A dehned by the hrst r rows (resp. last n — r rows). 
The boldface subscripts in A^ and A 2 are meant to distinguish these submatrices from the 
hrst and second row of A, to be written as Ai, A 2 . With this notation, the network dynamics 
is modelled by the DAE 

y' = —AiA'^^x 
= A2A^X. 


0 
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(10a) 

(10b) 


Splitting X = {y, z) in the right-hand side, and denoting Ljj = for ease of notation, 

this DAE can be rewritten as 


y' = -hiiy - L 12 Z (11a) 

0 = ^2iy + lh22^- (lib) 


Note that this constrained system has the differential-algebraic form displayed in ([3]). 

Theorem [1] describes the dynamics of the DAE flTTl) : in particular, it shows that the 
scheme defined above redistributes the commodity in a way such that all type-1 agents 
asymptotically accumulate the same amount of the resource. Note also that there is a 
somewhat unusual phenomenon, namely the existence of a line of non-isolated equilibria. We 
use from DAE theory the notion of a consistent initial value, which is a value of ( 2 /( 0 ), 2^(0)) 
which satisfies the constraints (lllbll (cf. [T ^ 130 1 Ii5]). 

Theorem 1. In the setting described above, the following assertions hold for the DAE / flil) . 

(a) The constraint llllb\] specifies an r-dimensional solution space. Its (transversal) inter¬ 
section with the hyperplanes Vi = ^ eM. defines a foliation of the dynamics by a family 
of (r — 1)-dimensional invariant spaces. 

(b) All eguilibrium points are located in the line yi = y 2 = ■ ■. = yr = zi = ... = Zn-r, which 
is transversal to the aforementioned invariant spaces if r >2. 


(c) The trajectory {y(t),z{t)) emanating from a given consistent initial value (|/(0),z(0)) 
converges exponentially to the eguilibrium 


Vi = Zj = 


2/1 (0) |/r(0) 

•) 

r 


1 , ■ ■ ■, A j = l,...n-r. 


Proof. 

(a) A well-known property in digraph theory states that any proper subset of the rows of the 
incidence matrix A of a connected digraph has maximal rank. This implies that the matrix 
A 2 , dehned by the last n — r rows of A, has maximal row rank (note that r > 1). In turn, 
this means that L 22 = A 2 AJ is a non-singular matrix; indeed, provided that A 2 A(^v = 0 we 
get v'^A 2 A^v = 0, and then Af^v = 0 yields n = 0. Since A 2 A 2 has order n — r, it follows 
that flllbp specifies an r-dimensional linear space. 

This space is hlled by solutions of the DAE, which are defined from those of the explicit 
ODE 


y' — — (Lii — Li2L22^1L<2i)|/ 


by means of the additional relation 

Z —— 11^22 lh2l^. 


( 12 ) 


(13) 
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That flllb|) intersects transversally the hyperplanes j/* = fc G M can be seen again as 
a consequence of the fact that L 22 = ^ 2^2 is non-singular, since the coefficient matrix of 
the set of linear equations dehning the intersection reads as 


1 , 


0 


IL21 IL' 


22 


where ij G has all entries equal to one. This matrix is easily seen to have maximal 
row rank, meaning that the intersection is indeed transversal and hence dehnes an (r — 1)- 
dimensional space, for any hxed k. 

Denoting by 1 the vector of I’s in R"^^, the invariance of these spaces follows from the 
identity A = 0 (expressing the fact that the sum of all rows of the incidence matrix A 
vanishes), which readily yields l'''L = 0 since L = AAA. Recasting flllbj) as —Lx = 0, from 
m one gets ijy' = 0, meaning that yi +.. .+yr indeed remains constant along trajectories. 


(b) The equations yi = y 2 = ■ ■ ■ = yr = A = ■ ■ ■ = A-r dehne a line of equilibria because 
of the identity = 0, which means that 1 spans kerL = AAA. Therefore any vector 
belonging to span < 1 > annihilates the right-hand side of fITT]) and therefore dehnes an 
equilibrium point of the DAE. There are no other equilibria because rk L = n — 1. 

Note that if r = 1, each one of the invariant spaces referred to in (a) amounts to a point, 
which actually belongs to the equilibrium line yi = Zi = Z 2 = ■ ■ ■ = z^-i- In cases with 
r > 2, the transversality of the invariant spaces mentioned above and the equilibrium line 
follows easily from the fact that, in R"", the hyperplane ~ ^ intersects transversally 

the line yi = ... = yr = Zi = ... = Zn-r- For later use it is worth emphasizing that the 
transversality of the intersection of the invariant spaces and the equilibrium line implies that 
each equilibrium is unique within each one of such invariant spaces. 


(c) The dynamics can be examined in R'" in terms of the ?/-variables via (IT^ . having in mind 
that the 2 ;-components are given by ([13]). In particular, the invariant spaces and the line 
of equilibria discussed in (a) and (b) are projected respectively onto the hyperplanes of R^ 
dehned by 

r 

^yi = k (14) 

i=l 


and the line span <lr>, that is 


yi = y2 = ■ ■ ■ = yr- (15) 

Notice that both spaces are orthogonal to each other. 

The rest of the proof relies on the fact that —(Ln — Li 2 L 22 ^L 2 i) is a symmetric, negative 
semidehnite, corank-one matrix. Its symmetry follows from that of Ln and L 22 , together 
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with the identity L 21 = 1^2 (recall that L^- = AiAj). To check that it is negative semidehnite, 
write 

-„yL„ - = -{v- - (J;;; J;-) ( ) < 0 

for any v & W, since —L = —AA^ is itself negative semidehnite. Similarly, the fact that 
— (Lii — Li 2 L 22 ^L 2 i) has corank one follows easily from the remark that 


V G ker (Lii — Li 2 L 22 ^IL 2 i) 


V 


IT . 


together with the identity corkL = 1. 

Altogether, this implies that —(Ln — Li 2 L^ 2 ^L 2 i) has a simple zero eigenvalne, its associ¬ 
ated eigenspace being the line of eqnilibria of flT^ . and r —1 real and negative eigenvalues; the 
corresponding eigenvectors are orthogonal to the equilibrium line because — (Ln—Li 2 L^ 2 ^L 2 i) 
is symmetric. Since the hyperplanes dehned by flMl) are themselves orthogonal to the equilib¬ 
rium line, the evolution in such invariant hyperplanes is characterized by these r — 1 negative 
eigenvalues. 

This means that trajectories evolve exponentially towards the unique equilibrium in each 
invariant hyperplane; it is a trivial matter to check that this equilibrium must be defined by 

J/l(0) HriO) . ^ 

Vi = -, t = 

r 

which is the unique solution to (lT4ll - (IT5|l with k = yi{0) + ... + |/r(0). Finally, the identities 


= 


1/1 (0) ?/r(0) 


J = 1, 


,n — r, 


follow from the fact that the equilibrium line spanned by in R*” is the ^/-projection of the 
line spanned by 1 in R”. □ 


2.3 Constrained dynamics on general networks 

The systematic analysis of constrained dynamics on general networks dehnes a topic of 
potential interest in many application fields. Just to name a few, the semiexplicit DAE (13|) 
may, in neural learning processes, model an equilibrium assumption for the dynamics of the 
neural state variables ( 2 ; in ([3])), whereas the evolution of the synaptic weights (corresponding 
to the //-variables in ([3])) characterizes the learning scheme; cf. [5[ll6]. Beyond ([3]), constraints 
may arise not involving the algebraic variables z] this yields a so-called Hessenberg DAE 

y' = f{y,z) (16a) 

0 = g{y). (16b) 


Such DAEs may also arise naturally in the presence of constraints involving only node 
variables (standing for y in flTHll ) in models which also use certain variables supported on 




the network edges (z in ffT6|l h Along the lines later discnssed in Section |3l Hessenberg 
DAEs may natnrally arise also in mnltilevel networks in which certain restrictions among 
gronp variables are not visible at the agent’s level, and viceversa. On the other hand, as 
it happens in nonlinear circnit theory (cf. [311 H5]). model redaction techniqnes in general 
networks may also lead to quasilinear DAEs of the form A(x)x' = f{x), where A{x) is an 
everywhere singular matrix-valued map. A reduction to an explicit ODE formulation may 
be difficult or impossible in these broader contexts. Besides circuit theory, where DAEs are 
now pervasive, different helds are starting to beneht from the systematic use of constrained 
formulations to describe network dynamics; examples are flow networks [27] (including flow 
dynamics in tree networks 1371), integrated process networks [15], dynamic traffic networks 
[21] , pressurized water networks [H] or, from a different perspective, multiagent descriptor 
systems [51]. Analytical results directed to general networks are very promising in these and 
other related areas. 

3 Dynamics on multilevel networks 

We now turn our attention to a dynamical problem in a clustered network. As detailed below, 
we will assume that the set of agents (nodes) is partitioned into disjoint clusters or groups; 
state variables will be dehned both on individual agents and on groups, and membership of a 
given group will entail a dynamical relation between the node and the group state variables. 
Additionally, we will assume that there are not only relations among individual agents, but 
also among groups. We are not concerned about the nature of these relations or the criterion 
according to which the nodes are clustered, but if it is of help the reader may think e.g. of 
a set of economic agents clustered in cooperatives, commercial transactions taking place 
both among agents and among cooperatives. Our goal is to analyze the dynamics of such a 
two-level network. 

3.1 Two-level networks 

We consider a two-level network dehned by a 4-tuple (V, E, H, Eh), where 

o 1/ is a set of n > 1 nodes or agents] 

o E is a set of m edges connecting some pairs of nodes; 

o H IS a family of p > 1 non-empty sets or groups of nodes; and 

o Eh is a set of q generalized edges connecting some pairs of groups. 

Note that both (V, E) and {H, Eh) are graphs, whereas (V, H) dehnes a hypergraph (with 
hyperedges corresponding to the above-dehned groups) if all nodes belong to at least one 
group [8]. This structure simply models a set of agents a) with a dyadic (binary) relation 
among them; b) joined together into certain groups; and c) displaying also a group-level 
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dyadic relation. Both (V, E) and {H, Eh) will be assumed to be simple graphs, and we will 


focus on cases in which H dehnes a partition of the set of agents into p pairwise-disjoint 


groups, the j-th one including nj > 1 agents. Without loss of generality we assume that the 
agents are numbered according to this partition, so that the indices 1 ,... ,ni correspond to 
the hrst group, rii -|- 1 ,..., ni -|- 77-2 to the second one, etc. 

We assume that a reference direction is given to each edge and each generalized edge, 
giving both (V, E) and (if. Eh) a digraph structure. The incidence matrix describing the 
dyadic relation at the agents’ level will be denoted by A, as in Section[2](cf. (jl])). Analogously, 
the incidence matrix describing the relation among groups will be written as Aq = (ap) G 

IRPXg, 



1 if generalized edge j leaves group i 


(17) 


—1 if generalized edge j enters group i 

0 if generalized edge j is not incident with group i. 


In turn, the entries of the incidence matrix Ah = (op) G of the hypergraph {V,H) 

read as 



1 if the 7-th agent belongs to the j-th group 
0 otherwise. 


(18) 


3.2 Node-group dynamics 

In the setting considered above, dynamics may take place both at the agent and at the group 
level; therefore, state variables will be dehned not only for agents but also for groups. For 
simplicity, we will assume that each agent and each group have exactly one (scalar) state 
variable, to be denoted by a:* for i = 1 ,..., n and pj for j = 1,... ,p, respectively. Alone, the 
X- and the 7/-variables might be understood to correspond to dynamical processes dehned on 
the graphs (y,E) and {H,Eh)- Certainly, the interest is placed on cases in which there is 
additionally a set of dynamic processes relating the agent state variables Xi with the group 
ones i/j. 

In this context, we wish to analyze an extension of the redistribution scheme discussed 
in Section [2] to two-level networks. Specihcally, both at the agent and at the group level 
there will be a dynamical process dehned by the corresponding Laplacian matrices, that 
is, L = AA^ and G = AqA^. Without an agent-group dynamic coupling, such dynamical 
systems would simply read as 


x’ = —Lx 


(19) 


and 


y' = -Gy, 


( 20 ) 


as in (E]). 
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The coupling between agents and groups will be defined by the assumption that the 
group variable yj stands for a collectively stored commodity at the j-th group, which (in the 
absence of transactions among agents or groups) evolves according to the individual amounts 
Xi stored by the agents which are clustered in that group. Specifically, we will assume that 
the rate at which the Tth agent contributes to yj (or collects from y^) is proportional to the 
difference 



The idea here is that the i-th agent contributes to (or collects from) yj depending on the 
difference between its own amount of the resource, Xj, and the part yj/uj that would corre¬ 
spond to each agent in an eventual uniform distribution of the group stored commodity yj 
to all rtj agents. This yields 


V • 

x'i = —Xi ^—for i in group j 

Hi 


and 

y'j = ^Xi-yj, j = l,...,p. 

i in j 

The latter relation is derived immediately from the identity 




i in j 


where for the moment we are disregarding the flows among agents and among groups. 

Note that, by construction, AJ^Ah is a diagonal matrix of order p, the j-th diagonal 
entry being Uj, that is, the number of agents in the j-th group. The node-group dynamic 
relations described above can be then globally expressed as 


x' = -x + Ah{A]jAh) ^y (21a) 

y' = ( 21 b) 


which can be understood as a dynamical process in the hypergraph (V, H). 

The whole dynamical process combines the flows among agents and among groups, de¬ 
scribed by flT^ and fl20|) . respectively, with the node-group dynamic relations defined by 
(1^ . This leads to the model 

x' = -(In + ^x + AH{A]jAH)~^y (22a) 

y' = A]jx - (Ip -F G)y, (22b) 


where and Ip are identity matrices of orders n and p. 
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3.3 A stability problem 


System fl2^ describes a dynamical process in a two-level network, withont any restriction so 
far in the topology of the graphs (V, E) and (if, Eh)- We address here a stability problem in 
a simplihed setting; we assnme that there are no interactions at the agents’ level or, in graph- 
theoretic terms, that all nodes are disconnected. The flow (l20ll describing the interactions 
among gronps in the graph [H, Eh) and the node-gronp dynamics fl2T]l in the hypergraph 
(1/, H) still apply, to yield the model: 

x' = -X + AH{A]jAH)~^y (23a) 

y' = - (I + G)y. (23b) 

Below we denote by x(i) the gronp to which agent i belongs. 

Theorem 2. Consider system (d^) and assume that the graph {H,Eh), describing the rela¬ 
tions among groups, is connected. 

(a) System has a line of eguilihria defined by the relations 

Xi = (i = 1,..., n), = 7/2 = ... = 2/p. (24) 

n-xii) 

(b) An initial condition (a:(0),//(O)) converges exponentially to the eguilibrium defined by 

1 / n p \ 

Xi = i = l,...,n, yj = Y I + ) ’ J = (25) 

P \i=i k=i / 

Proof. 

(a) By means of a Schnr redaction of fl23|) . eqnilibria are easily checked to satisfy 

A]jAh{A]jAh) ^y — (IG)y = 0, 


that is. 


Gy = 0. 


Since the graph {H,Eh) is connected, G = AqA^ is rank-deficient by one, and kerG = 
ker Aj is dehned by the relations 2/i = 2/2 = • • • = 2/p- The x-components are dehned by 
X = AH{AljAH)~^y. By constrnction AffiAn is a diagonal matrix, and the j-th diagonal 
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equals rij, that is, the number of agents in the j-th group; we then have 


Ah = 


/ 

1 

0 

0 

0 

... 0 \ 





0 

0 

0 

... 0 




(’^l) 

0 

0 

0 

... 0 





0 

0 

0 

... 0 


( in 

0 


1 

0 

0 

0 

... 0 


ni 

0 

J_ 


0 

1 

0 

0 

... 0 


n 2 


0 


0 

0 
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( 26 ) 


so that X = AH{AljAH)~^y amounts to Xi = for i = 1,..., n, as stated in 

(b) System will be proved stable by using Gersgorin theorem [25], according to which 
all eigenvalues of the coefficient matrix of the right-hand side of fl22]) . namely 


M = 


-I Ah{AJjAh) ^ 
-(I + G) 

must lie on the union of the discs \z — frijjl < Rj, with 


4 T 


n+p 






i=l 


We are denoting by rriij the entries of the matrix M. It is obvious that 


rrijj = —1 if 1 < j < n, 


and 

TJijj = —1 — dj-n ifn + l<j<n + p, 

where dj stands for the number of connections of the j-th group to other groups (that is, its 
degree in the graph {H, Eh))] we are making use of the fact that G is the Laplacian matrix 
of the graph {H, Eh) describing the group connections. 

On the other hand, there is only one non-vanishing entry (with value 1) in any column 
of AJhi so that 

= 1 if 1 < j < n. 
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This means that n discs are centered at —1 and have radius 1. Additionally, by construction 
one can see that the sum of the entries in any column of is 1; indeed, from 

fl26|l we may check that the j-th column of Ah{AA^Ah)~^ has rij nonvanishing entries equal 
to l/rij. Together with the fact that the off-diagonal entries of —G define the adjacency 
matrix of the graph (if, G^), this yields 

Rj = 1 + dj-n foTn + l<j<n + p 

so that the remaining p discs are centered at —1 — dj^n and have radius 1 -|- dj^n (with 
j = n + 1,... ,n + p). For later use, notice that this also shows that the sum of the entries 
of any column of M does vanish. 

The remarks above prove that all Gersgorin discs are located on the left half complex 
plane, except for a tangency with the imaginary axis at the origin. From (a) it follows that 
the matrix M has a unique zero eigenvalue (with eigenvectors in kerG), so that the other 
n + p — 1 ones are away from the imaginary axis. 

The last remark needed to prove (b) is that the quantity Vj preserved 

along trajectories, which is a consequence of the aforementioned fact that the sum of the 
entries in each column of M is null (so that + ^i=i d'j ~ ^)- This means that 

n p n p 

^ Xi + ^ Xi(0) + vM 

i=l j=l i=l j=l 

is a family of invariant hyperplanes, the evolution rate in each one of them being characterized 
by the n-|-p —1 eigenvalues which are located on the left half-plane. Therefore, the dynamics 
evolves exponentially towards the point defining the intersection of such invariant planes and 
the line of equilibria arising in (a). By combining the identities yi = y 2 = ■ ■ ■ = Up with the 
fact that the conditions Xi = y^(i) / n^{i) yield Xi = Y7j=i Vj equilibrium, it is easy to 
check that this intersection is given by fl25|l and the proof is complete. 

□ 


Theorem [2] extends the results discussed in Section [2] to the two-level dynamic network here 
considered. Now the scheme converges asymptotically to an equal distribution of the resource 
among all groups, and to a uniform distribution among the individual agents in each group. 
Differences may arise between agents of different groups. 

We leave for future study the full analysis of (|22D . which may display more intricate 
dynamics; further research may extend these results to problems with flowrates different 
from the ones yielding the graph Laplacian matrices in flT^ and (|20|) . or with node-group 
dynamic relations different from Our results should also be of interest in an eventual 
deterministic, continuous-time dynamic analysis of multilevel networks in general, possibly 
accommodating more aggregation levels. 
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4 On the equilibrium set of potential-driven flow networks 

Theorems [U and [2] above share a somewhat unusual property in dynamical systems theory, 
namely the existence of a line of equilibria. Non-isolated equilibria have received attention 
within the theory of bifurcation without parameters (cf. [33] and references therein; see 
also |1]). In our context, the key element supporting the existence of a line of equilibria 
is the fact that systems flTTl) and fl2^ essentially describe the dynamics of potential-driven 
flows. As detailed later, equilibria dehne a one-dimensional manifold (a line) because of 
the implicit assumption that the networks involved are positively weighted. Without this 
assumption, the structure of the equilibrium set may be more intricate; in this section we 
examine such structure in flow networks without this positiveness assumption, arriving at 
a graph-theoretic characterization of the structure of the equilibrium set in terms of the 
network spanning trees. In particular, we will apply our results to networks based on signed 
graphs, originally introduced by Harary izg. 

4.1 Nonlinear flows 

A potential-driven flow on a (directed) graph with n nodes and m edges is deflned by a 
dynamical system of the form 


aix[ = -Aif{A^x), i = 1,..., n. (27) 

Here x = (xi,..., x„) and Ai is the Ath row of the incidence matrix A as deflned in (jlj). In 
turn / : —)■ is a possibly nonlinear, differentiable map describing the flowrates in the 

edges. It is assumed to have a diagonal structure, that is 

/ = /i X ... X /^, 

with /j : M —)■ M depending only on {A'^)jX, where {A'^)j is the j-th row of A'^] this means 
that the flowrate in edge j (connecting nodes i and k) depends only on the difference Xi — Xk- 
The variables x can be thought of as a potential and hence the “potential-driven” label for 
the flow: electric potential or pressure are examples in electrical circuits and water networks, 
respectively. Finally, the coefficients a* are 0 or 1 for each node; the case «* = 1 (resp. = 0) 
corresponds to nodes which accumulate (resp. do not accumulate) a certain amount of the 
quantity or resource (e.g. electrical charge, water, gas, a given commodity, etc.) which is 
flowing in the network; note that 1^7} describes the continuity equations at both types of 
nodes. 

A simple example with a, = 1 for all nodes, / being the identity map, is deflned by 
(IS])-(l6l), yielding the Laplacian dynamics ([7]). A problem combining nodes with = 1 and 
ttj = 0 is given by system CD. analyzed in Section |2l Additionally, the reader may think of 
a (possibly nonlinear) resistive circuit as an example with = 0 for all nodes. 
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4.2 Equilibria and the subimmersion theorem 

To make things simpler, in the sequel we assume that Oj = 1 for all nodes, focusing on 
dynamical systems of the form 

x' = —Af{A^x), (28) 

the aforementioned restrictions on the form of / still holding. Also for the sake of simplicity 
we assume that the digraph is connected, so that ik A = n — 1. For notational convenience 
we denote by F{x) the right-hand side of (1251) . that is, 

F{x) = -Af{A^x). (29) 

It is easy to check that equilibria of fl25D may never be isolated; indeed, provided that 
F{x) vanishes for a given x*, then x = x* + v also annihilates F{x) = —Af{A'^x), for any 
V G kerA'''. Note that this kernel is never trivial, being one-dimensional in a connected 
digraph. The problem we address in this section is the characterization of the cases in which 
the equilibrium set is locally a line, as it happens in Theorems [T] and [2J We show below that 
the one-dimensional nature of the equilibrium set in these theorems is implicitly supported 
on a positiveness assumption on the digraph weights, without which the equilibrium set may 
locally have a higher dimension. In general, we do not require / to be linear; when / is a 
linear map the equilibrium set is a linear manifold and the results hold globally. 

We will make use of the subimmersion theorem, which states that if hi is an open subset of 
M"’, and F is a smooth mapping hi —)■ such that the Jacobian matrix F'{x) has constant 
rank r < p on Q, then for every y in F{Vt) the set F~^{y) is a submanifold of hi with 
dimension n — r (see e.g. Th. 3.5.17 in pQ or Th. III.5.8 in HU]). The result also holds if F 
has constant rank r on a neighborhood of F~^{y). We will use this result with F given in 
(EH), y = 0, p = n and r = u — 1 to characterize the problems in which the equilibrium set 
is locally a line around a given equilibrium point x*. The Jacobian matrix F'{x) reads as 

F'{x) = -Af\A'^x)A'^ (30) 

and, since rk A = n — 1, it follows that F'{x) is persistently rank dehcient. Therefore, for 
the equilibrium set of flUU]) to be one-dimensional (at least locally around x*), it is enough to 
derive conditions guaranteeing rk F'{x*) = n — 1, since this maximum rank would necessarily 
be attained also on a neighborhood of x*. In order to examine the rank of F'{x*), let us 
denote the derivatives of the components of / at x* as 

W, =/'(AV), J = l,...,m, (31) 

and let W stand for the diagonal matrix with entry Wj in the j-th diagonal position. Note 
that this matrix is indeed diagonal because of the assumption that fj depends only on 
{A'^)jx. With this notation we have 

F'{x*) = -AIFA'^. (32) 

This expression shows that the Jacobian matrix F'{x*) has the structure of a weighted 
Laplacian matrix, with weights being defined by the derivatives fj{A'^x*). 
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4.3 Positive weights 

If all weights Wj (that is, all derivatives fj{A^x*)) are positive, then it is a simple matter to 
check that 

kei AWA^ = kei . (33) 

Indeed, just note that AWA'^u = 0 implies u^AWA'^u = 0 and therefore AAu = 0 because 
of the positiveness of the diagonal matrix W. The relation depicted in fl33p implies that 

rk F'{x*) = rk AWA'^ = rk A'^ = n - 1 

and, because of the subimmersion theorem, it follows that equilibria actually dehne a curve 
near x*. 

Implicitly, this underlies the existence of a line of equilibria in the setting of Theorems 
[T] and [21 Note that system ([7]), dehned by the Laplacian matrix AA^, can be understood 
as a particular instance of the product AWA^ with W = Im- Actually, one may show that 
the results of Sections [2] and [3] actually hold if the product AA^ is replaced by AWA^ with 
a positive, diagonal W] in other words, those results still apply if the networks are assumed 
to be positively weighted. Note that the linear setting considered in those sections yields a 
linear manifold of equilibria and avoids the need for a local approach. 

4.4 Negative weights and the structure of the equilibrium set 

If W in (152]) includes negative entries, the remarks just stated do not apply, and we may 
actually hud potential-driven flow dynamics on connected graphs displaying higher (> 2) 
dimensional manifolds of equilibria. This may be the case even in linear problems: examples 
in a linear context, involving signed graphs, can be found in subsection 14.51 below. 

A different approach is needed to analyze the structure of the equilibrium set in digraphs 
with possibly negative weights. In this setting, we will drive the problem to a context known 
in the framework of circuit theory, namely the analysis of nodal admittance matrices, which 
can be performed along the topological approach stemming from the work of J. C. Maxwell 
[m 031 . In Theorem 121 which is the main result of this section, we make use of the notion 
of the weight of a spanning tree; this is simply the product of the weights of the tree edges. 

Theorem 3. Assume that the dynamical system is defined on a connected digraph. The 
Jacobian matrix F\x*) in / f72l) has corank one if and only if the sum of the spanning tree 
weights does not vanish, where the edge weights are given by / fTTl) . If this sum is not null, 
then the set of eguilibria is a curve locally around the eguilibrium point x*. 

Proof. The matrix AWA^ in the right-hand side of (1321) is a weighted Laplacian matrix, 
which is known to be rank dehcient since rk A = n — 1. The hrst step in the proof shows 
that all (n — l)-minors of AWA^ do vanish if and only if a single one of them does; therefore, 
the identity rk AWA'^ = n — 1 will hold if and only if one (hence any) principal minor does 
not vanish. 
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By construction, minors of AWA'^ of order n — 1 are determinants of products of the form 

Ar.WAl, 

where A^^ and A^^ are reduced incidence matrices, dehned by two arbitrary choices of {n — 1) 
rows of A. The key remark here comes from a graph-theoretic property, saying that any 
set of n — 1 rows of the incidence matrix of a connected digraph are linearly independent 
(actually dehning a basis of the so-called cut space] cf. [S]). This means that a relation of 
the form 

A,^ = KAr, 

holds for a non-singular matrix K. Therefore 

det{A,^WAjj = det(klrilTklJj det itT 

and it follows that detAjj does vanish if and only if det(y4rihTACj does. Hence, the 
eventual vanishing of all (n — l)-minors occurs if and only if a single one of them is null; 
equivalently, in order to examine the condition rk AWA^ = n — 1 it suffices to study the 
vanishing of a principal minor, that is, the condition 

det{A,WAj) ^0, (34) 

for an arbitrary choice of a reduced incidence matrix A,.. Without loss of generality, in what 
follows we are allowed to work e.g. with the reduced incidence matrix A,, dehned by the 
hrst n — 1 rows of A. The corank-one condition on F'{x*) then amounts to evaluating the 
condition fl54D . 

This can be performed using a Cauchy-Binet determinantal expansion of fl34p (cf. |26j). 
which yields 

dei{A,WAl) = ^det(H")det(W^)det((H7)^), (35) 

the sum being taken over all index sets a, fd F {1 ,... ,m} with n — 1 elements. In the 
submatrices involved in this expansion, the subscripts (resp. superscripts) a, {3 are used to 
select a set of rows (resp. columns); for instance, denotes the submatrix of specihed 
by all rows and the columns indexed by a, whereas is the submatrix of W dehned by 
the rows and columns specihed by a and /5, respectively. 

The expansion depicted in 0351) can be simplihed using the following two remarks. First, 
a set of columns of the reduced incidence matrix is known to yield a non-singular matrix 
if and only if such columns correspond to the edges of a spanning tree; in that case, we have 
det(y4") = ±1, because incidence matrices are totally unimodular (see e.g. [3]). 

Second, because of the diagonal form of hF, one can easily check that det{W^) does 
vanish ii a ^ [3] when a = /3, we have 

det(W'“) = n 

jea 
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which is the tree weight when a specihes a spanning tree. Note additionally that, li a = (3 
then det(A") = ±1 eqnals det((A7)/3) = clet((A7)a)) so that 

det(74") det((74j'')o) = 1. 

Altogether, these remarks prove that the determinantal expansion displayed in (135|1 
amounts to 

det (ArWAj) = 

a jGa 

where the sum is taken over the index sets a which specify a spanning tree. It follows that 
any (n — l)-minor of the Jacobian matrix F'{x*) is non-null, and therefore F'{x*) has corank 
one, if and only if the sum of the spanning tree weights does not vanish, as claimed. 

Provided that the non-vanishing condition flM|) holds, the identity rk F'{x*) = n — 1 
yields, as a direct consequence of the subimmersion theorem in the terms stated above, a 
local one-dimensional structure for the equilibrium set near x*, and the proof is complete. 

□ 


Theorem [3] provides a graph-theoretic characterization, in terms of the digraph tree struc¬ 
ture, of the problems which systematically lead to curves (and not higher dimensional or 
singular manifolds) of equilibria in potential-driven flow networks with (possibly) some neg¬ 
ative weights. Stability aspects in this context are in the scope of future research; note in 
particular that the failing of the non-vanishing requirement in the sum of tree weights might 
be responsible for bifurcation phenomena in nonlinear problems. 

Certainly, if all weights are positive then the sum arising in Theorem [3] is positive and 
therefore non-null, because all tree weights are positive. This is of course consistent with 
the discussion of subsection 14.31 regarding the one-dimensional nature of the equilibrium set 
in positively weighted networks, for which there is no need to use these tree-based tools. 

4.5 Signed graphs 

A nice corollary of Theorem |3] holds for signed graphs, originally introduced by Harary [22] 
and widely used in balance theory (cf. [23] and more recent references such as [HI [I6]). A 
signed graph or s-graph is a graph iV,E) endowed with a map E —)■ {—1,1}, that is, an 
assignment of either a -|-1 or a —1 weight to all m edges. As above, we will let W G 
stand for the diagonal matrix of weights. 

In this context, a flow problem which arises as a particular case of fl28D is 

x' = -AWA'^x. (36) 

This system can be again understood as a flow, via the continuity equations x' = —Au, 
in which the flowrates u are given by the linear relation WA'^'x. With respect to (ISil-lp]), 
yielding the Laplacian dynamics (|7]), now the presence of —1 weight values models a flow in 
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which adjacent agents i, k tend to increase the difference between Xi and Xfc, since the ffow 
from i to k now equals Xk — Xi (instead of Xi — Xk, as in subsection 12 .ip . 

Contrary to the results in Sections [2] and [3|, in which equilibria define a line, the presence 
of negative weights might now result in higher dimensional equilibrium manifolds. The cases 
in which this may happen are exactly characterized in Corollary [1] below. Since all weights 
are +1 or —1, we may now define a spanning tree as positive or negative simply if its weight 
product is +1 or —1, respectively, or, equivalently, if it contains an even (resp. odd) number 
of edges with negative weight. 

Corollary 1. Let the dynamical system defined on a connected signed digraph. Then, 

the dimension of the eguilibrium set is higher than one if and only if the numbers of positive 
and negative spanning trees coincide. 

This results follows immediately from Theorem [3] since, for the sum of weight products 
to vanish in a signed graph, the amount of positive trees (which are responsible for a +1 
term in the sum) must obviously match the number of negative trees (which yield a —1 term 
in the sum). 


Examples. Simple examples illustrating the result above can be defined using the graphs 
shown in Figure [TJ The first case, on the left of the figure, is simply a 4-cycle with two 
positive signs (continuous lines) and two negative signs (dashed lines). This graph has just 
four spanning trees, two of which are positive and the other two negative. The second 
example, on the right, depicts a complete graph with two positive and four negative 
signs; this is a so-called balanced structure [22], in which the nodes may be split in two 
groups in a way such that all edges inside a group have positive weights and all connections 
between groups are negative. In this example the groups are defined by the two nodes on 
the left and those on the right, respectively. According to Cayley’s formula this graph has 
4^ = 16 spanning trees, and it is not difficult to check that exactly half of them are positive. 
This means that the dynamical system (1361) should exhibit in both cases a linear manifold 
of equilibria with dimension greater than one. 

In the first case, we number and orientate each edge beginning on the top and according 
to a clockwise orientation of the cycle. This yields 

0 0 -1 \ 

1 0 0 

-1 1 0 ■ 

0-1 1 y 


/ 1 
-1 
0 

V 0 


With W = diag (1,1, —1, —1), the right-hand side of (1361) reads in this case as 


-AWA'^ 


/ 0 1 0 -1 \ 
1-210 
0 10-1 

^ -1 0-1 2 y 
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Figure 1: Sign assignments yielding degenerate flow dynamics on (a) a 4-cycle; (b) K^. Edges 
with a negative sign are dashed. 


The kernel of this matrix defines the equilibrium set and is defined by the relations 

a:i + X 3 

X2 = Xi = - - -, 

hence defining a two-dimensional linear manifold (a plane) of equilibria. Notice that these 
equilibrium points arise from a constant flow Xj — Xj+i (i = 1,..., 4, with the terminological 
abuse X 5 = Xi) which annihilates all derivatives x', flowing in the clockwise (resp. counter¬ 
clockwise) direction if xi > X 3 (resp. if Xi < X 3 ). Note that the equilibrium plane includes 
the line Xi = X 2 = X 3 = X 4 for which the flow vanishes. 

In the example defined on (on the right of Figure [1]), we number and orientate the 
edges according to the following incidence matrix: 

( I 0 0-1 1 o\ 

-1 1 0 0 0 1 

0 - 110 - 10 ’ 

^0 0-1 1 0 -1 y 

and the weights W = diag (—1,1, —1,1, —1, —1) yield 

1-1 l\ 

1 1-1 
1 1 -1 ■ 

1-1 ij 

Now the equilibrium set is three-dimensional, being defined by the identity 


-AWA^ = 


( 1 - 
-1 
-1 

V 1 - 


xi -F X 4 = X 2 + X 3 . (37) 

Again, these equilibrium solutions yield nonvanishing flows in the graph edges, except for 
those in the line xi = X 2 = X 3 = X 4 . This example also illustrates that the rank drop in the 
product AWA^ may be greater than two. 
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We leave for future work the analysis of the nature of such equilibrium solutions exhibiting 
non-vanishing flows, that is, equilibria which are not in the line xi = X 2 = ... = Xn] this 
analysis should explain, for instance, the relation of the equilibrium solutions in the second 
example with the presence of a balanced structure in the graph: note that flTTIl expresses that 
the total amount of the resource x stored in both groups (namely, {xi^x^} and {x 2 ,x^}) is 
the same. More generally, further research might address stability aspects and bifurcations in 
these negatively weighted networks, not only in signed graphs but specially in the nonlinear 
setting introduced in subsection 14.11 

Note hnally that even though in Section 0] we have focused for simplicity on flow networks 
with homogeneous agents (because of the assumption ctj = 1 for all nodes, cf. subsection 
14.2p . the results are also of potential interest in broader contexts, involving e.g. heterogeneous 
nodes or group dynamics as in Sections [2] and [3l 

5 Concluding remarks 

We have addressed in this paper some qualitative problems involving network dynamics be¬ 
yond the scenario dehned by a set of dynamical systems supported on the nodes of a graph. 
First, the attention has been focused on potential-driven flow dynamics either on networks 
with heterogeneous agents or on multilevel networks. Differential-algebraic models, arising 
here in flow networks with heterogeneous agents, are worth receiving further attention in or¬ 
der to model and analyze constrained dynamics on general networks. Dynamical systems on 
multilevel networks also have a broad scope for future research. From a different perspective, 
we have analyzed the structure of the equilibrium set in (possibly nonlinear) flow networks 
with negative weights. The results in this regard apply in particular to signed graphs, and 
the connection of our results to the theory of balance and clusterability in signed graphs re¬ 
quires further analysis. In greater generality, stability properties and bifurcations are worth 
being studied in networks with negative weights. 

Our research is of potential interest in the analysis of social and economic networks 
involving flows. Future work might also extend these results to problems displaying network 
evolution, and also to stochastic and/or discrete-time contexts. 
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